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I.  INTRODUCTION 


The  objective  of  the  work  described  in  this  report  is  to  develop  a 
computer  program  for  modeling  the  seismic  event  identification  capability  of  an 
existing  global  network  of  seismograph  stations.  In  addition  the  program  is  to 
provide  for  evaluation  of  the  impact  on  the  overall  network  identification 
capability  resulting  from  modifications  to  the  network  in  the  form  of  additions, 
deletions  or  geographic  relocations  of  stations. 

Early  work  on  this  topic  was  conducted  at  the  Defense  Advanced 
Research  Projects  Agency  (DARPA)  by  Drs.  Ralph  W.  Alewme  and  Carl  F. 
Romney.  Based  on  the  results  of  their  efforts  it  was  decided  to  proceed  with  the 
design  and  implementation  of  an  interactive  program  that  would  achieve  the 
above  stated  objectives. 

The  initial  approach  adopted  was  to  construct  an  algorithm  that 
combines  a  statistical  analysis  of  existing  discrimination  data  with  analytical 
network  detection  and  location  modeling  techniques.  More  specifically, 
identification  probability  versus  m^  curves  for  an  existing  network  were  to  be 
estimated  from  an  analysis  of  past  discrimination  performance.  In  order  to 
predict  the  identification  curves  for  a  modified  network,  the  curves  for  the 
existing  network  had  to  be  changed  to  reflect  expected  differences  in  the 
detection  and  location  capabilities  of  the  existing  versus  the  modified  network. 

After  considerable  effort  the  analytical  approach  proved  to  be  intractable 
m  large  part  owing  to  the  complex  dependence  of  seismic  discriminants  on  the 
varied  properties  of  seismic  phases  used  in  the  application  of  these 
discriminants  as  well  as  the  number  and  distribution  of  detecting  stations 
Subsequently,  based  on  suggestions  by  Dr.  Robert  R  Blandford  of  DARPA.  a 
different  approach  (Monte  Carlo)  was  adopted  and  implemented  in  the  network 
identification  modeling  program  The  resulting  program  is  referred  to  as  NICE 
(Network  Identification  Capability  Evaluation) 
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In  the  Monte  Carlo  approach  the  source  properties  of  seismic  events  in 
a  geographic  region  of  interest  are  generated  by  random  selection.  The  NICE 
program  then  computes  synthetic  seismograms  of  these  events  for  a  specified 
network  of  seismograph  stations  and  determines  discrimination  scores  based 
on  the  synthetic  seismograms.  This  approach  has  several  strong  advantages, 
the  most  important  of  which  is  that  the  elements  of  the  problem  affecting 
network  identification  capability  are  included  naturally  and  directly  and  can  be 
tested  independently  for  accuracy.  The  elements  are  divided  into  three  major 
categories: 

1.  source  and  propagation, 

2.  network-controlled. 

3.  processing. 

The  methodology  incorporated  in  the  NICE  program  allows  for  the 
straightforward  inclusion  of  such  effects  as  the  reduction  of  surface  wave 
amplitudes  with  event  focal  depth  (and  the  corresponding  reduction  in  Ms:mb), 
the  shift  in  corner  frequency  with  source  dimension  (and  corresponding  changes 
m  Mg  and  m^),  the  relative  frequency  content  of  P  and  S  waves,  focal 
mechanisms  and  source  type  (explosion,  earthquake,  spall  model,  etc  ).  In 
addition,  because  of  the  modular  structure  of  the  program  new  identification 
criteria  can  be  readily  included  Should  the  new  criteria  require  that  seismic 
phases  not  currently  computed  in  the  synthetic  seismogram  section  of  the 
program  be  included,  they  may  also  be  added  in  a  straightforward  way  In 
general,  the  program  has  been  structured  for  use  in  an  interactive  mode  with  as 
much  flexibility  as  possible  The  NICE  program  is  operational  on  a  Gould  6031 E 
computer  at  S-CUBED  and  on  a  VAX  1 1/780  at  the  DARPA  Center  for  Seismic 
Studies  (CSS)  m  Arlington,  Virginia. 

The  following  section  of  this  report  gives  a  detailed  description  of  the 
structure  and  flow  of  operations  in  the  NICE  program.  In  the  final  section  a 
Users'  Manual  is  given  to  facilitate  operation  of  the  program 


II.  THE  NICE  PROGRAM 


The  NICE  program  consists  of  four  major  elements  which  are  executed 
in  series,  as  shown  in  the  flow  chart  in  Figure  2.1.  The  overall  approach  is  a 
Monte  Carlo  simulation  in  which  events  are  generated  from  random  series, 
‘‘recorded"  in  a  network  as  synthetic  seismograms  and  then  "analyzed"  for 
discrimination  scores.  The  first  element  is  the  module  which  generates  a  series 
of  events  whose  statistical  properties  match  those  of  observed  events.  The 
second  module  computes  amplitudes  and  travel  times  for  these  events  using 
synthetic  seismograms.  These  data  are  analyzed  in  the  third  module  and  given 
discrimination  scores.  The  output  of  this  module  is  a  bulletin  which  resembles  a 
discrimination  bulletin.  Finally,  the  performance  of  the  discriminants  as  applied 
to  the  synthetic  bulletin  can  be  examined  using  plots  and  lists  displayed  by  the 
fourth  module.  The  characteristics  of  the  synthetic  bulletin  can  also  be 
compared  with  actual  observations  of  the  other  networks  using  the  fourth 
module.  In  the  following,  we  briefly  describe  these  four  modules.  Details  of  the 
calculations  are  given  in  the  appendices. 

2.1  THE  EVENTS  MODULE 

The  events  generator  creates  a  list  of  earthquake  epicenters,  depths, 
origin  times,  scalar  moments,  and  fault  orientations.  The  values  in  the  list  are 
members  of  statistical  ensembles  whose  parameters  have  been  deduced  from 
observations.  We  have  endeavored  to  formulate  the  program  so  that  each  set  of 
parameters  can  be  checked  independently  against  an  observable.  We  have 
allowed  for  regionalization  of  the  parameters,  a  feature  that  is  very  important  for 
modeling  the  differences  in  the  performance  of  discriminants  in  different  tectonic 
regions  of  the  USSR.  For  example,  the  distribution  of  events  with  depth  in  the 
Kurile-Kamchatka  area  differs  greatly  from  that  in  the  Hindu  Kush,  and  thus  the 
depth  discriminants  behave  differently. 


STATION  FILES 


Figure  2.1.  Flow  chart  of  Network  Identification  Capability  Evaluation  (NICE)  program. 
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Epicenters  are  distributed  uniformly  within  boxes  on  the  surface  of  the 
globe.  The  sides  of  the  boxes  are  great  circles  intersecting  at  right  angles  The 
user  specifies  two  latitude,  longitude  pairs  and  the  length  of  the  side  normal  to 
the  line  connecting  the  pairs.  The  formulas  for  finding  the  epicentral  coordinates 
are  given  in  Appendix  A.  Maps  made  from  synthetic  distributions  of  events 
specified  in  this  way  are  very  similar  to  the  actual  distribution  of  events  reported 
by  the  Soviets  to  the  International  Seismological  Center.  We  find  that  five  boxes 
provide  a  good  descriotion.  with  the  boxes  encompassing  the  seismic  areas  of 
the  Caucasus,  the  Hindu  Kush,  far-eastern  Siberia,  Tibetan  border-Lake  Baikal 
area,  and  Kurile-Kamchatka. 

The  distribution  of  focal  depths  is  specified  by  a  piece-wise  linear 
function  which  is  the  cumulative  number  of  events  occurring  at  or  below  a  given 
depth.  Appendix  A  gives  the  algorithm  for  focal  depth.  The  distributions 
observed  on  the  Soviet  bulletin  are  fit  well  by  this  specification,  often  with  a 
single  linear  segment. 

Fault  orientation  is  described  by  the  strike,  dip  and  rake  of  the  slip 
vector.  Each  quantity  is  distributed  uniformly  within  a  standard  deviation  of  a 
mean  value  In  most  seismic  areas,  a  few  faulting  mechanisms  dominate, 
making  the  choice  of  parameters  apparent.  In  areas  where  there  are  extensive 
seismic  determinations  of  the  fault  plane  (such  as  done  by  Murphy  and  Lambert. 
1985,  for  the  Hindu  Kush  and  Pamirs),  we  have  reliable  values  for  the 
distribution. 

The  recurrence  interval  of  events  at  a  given  scalar  moment  is  assumed 
to  follow  a  Poisson  distribution,  as  suggested  by  Lomnitz  (1966)  The 
cumulative  number  of  events  occurring  per  year  at  or  above  a  given  magnitude 
is  predicted  by  this  distribution.  This  is  typically  described  for  magnitudes  with 
"a"  and  "b"  values  (intercept  and  slope  on  a  number  versus  magnitude  plot) 
Analogous  values  for  log  moment  are  used  by  EVENTS  to  determine  the  origin 
time  and  scalar  moment  of  each  event.  See  Appendix  A  for  details  As  it  would 
be  tedious  and  perhaps  error  prone  to  require  the  user  to  enter  these 
parameters  interactively  and  since  the  parameters  are  not  changed  often,  the 
EVENTS  program  reads  files  containing  the  parameters  specific  to  each  seismic 
area. 
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2  2  THE  SYNWAV  MODULE 

The  SYNWAV  program  reads  the  list  of  source  characteristics  created  by 
EVENTS,  computes  synthetic  body  and  surface  wave  seismograms  for  these 
events  at  the  stations  in  the  network,  and  measures  the  amplitudes  and  travel 
times  from  the  synthetic  records.  These  results  are  used  to  construct  a  station 
bulletin  which  will  be  used  in  the  next  module  in  the  series,  ANLYST.  Details  of 
the  formulation  are  given  in  Appendix  B. 

The  first  step  is  to  read  several  files  containing  data  used  in  the 
calculations.  This  includes  a  station  file  which  includes  station  locations  and 
long  and  short  period  noise  parameters.  Amplitude  and  travel  time  tables  are 
read  for  the  body  wave  calculations  and  dispersion  and  eigenfunction  files  are 
read  for  surface  waves. 

The  body  wave  calculations  begin  with  computing  a  time  series  which  is 
the  convolution  of  the  instrument  response  and  t*  attenuation  operator  for  each 
phase  (P,  pP,  S).  We  assume  at  this  stage  that  propagation  introduces  no  other 
pulse  distortion.  The  body  wave  calculations  then  follow  these  steps  for  each 
event. 

1.  The  azimuth  and  range  from  the  source  to  each  station  is 
computed. 

2.  The  amplitude  decay  rates  and  travel  times  are  looked  up 
in  the  tables  giving  the  values  as  functions  of  distance  and 
focal  depth.  The  ray  parameters  are  found  by  numerically 
differentiating  the  travel  time  curves  with  respect  to 
distance  (Equation  Bi.3). 

3.  A  frequency-independent  factor  due  to  each  ray's  position 
on  the  focal  sphere  is  computed  assuming  that  the  source 
radiates  from  a  point  (ppmt  double-couple  approximation). 

These  are  the  factors  aK.  etc.,  in  Equation  (Bi  4). 

4.  The  instrument  t*  time  series  is  convolved  with  a  moment 
function  in  Equation  (Bl.5)  The  convolution  is  scaled  by 
the  radiation  pattern  and  decay  rate 
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5.  The  peak-to-peak  amplitudes  and  periods  are  measured'' 
for  the  a,  b,  and  c  swings  on  the  record. 

6.  The  noise  preceding  each  phase  is  found  from  the 
maximum  peak-to-peak  value  in  a  random  series  of 
duration  10  seconds  with  mean  and  standard  deviation 
prescribed  by  the  station  file.  In  addition,  noise  with  the 
same  properties  is  added  to  the  amplitude  measurements. 

T  or  the  body  wave  calculations,  we  assume  that  the  major  part  of  the 
propagation  path,  the  mantle,  introduces  no  pulse  distortion  (the  mantle  transfer 
function  has  a  flat  spectrum).  For  surface  waves  however,  dispersion  causes  the 
shape  of  the  wave  tram  to  change  significantly  with  range.  Thus  we  must 
compute  a  seismogram  for  each  station-event  pair,  rather  than  scaling 
amplitudes  as  is  done  for  body  waves  The  most  time  consuming  part  of  a 
surface  wove  calculation  is  finding  the  dispersion  curves  and  corresponding 
eigenfunctions  for  the  fundamental  mode  To  reduce  computation  time,  the 
group  and  phase  velocities  for  the  chosen  propagation  path  are  computed  by 
the  SYNSRF  program  (a  general  purpose  surface  wave  program  written  at 
S-CUBED)  for  a  chosen  propagation  path  and  stored  on  a  file.  In  addition,  the 
eigenfunctions  for  a  suite  of  source  depths  are  calculated  and  stored  SYNWAV 
then  reads  the  eigenfunctions  and  interpolates  for  the  values  of  source  depth 
required  by  the  events  catalog.  The  interpolated  values  are  then  included  m  the 
surface  wave  formulation  shown  in  Appendix  B.  The  peak  value  whose  period  is 
between  18  and  22  seconds  is  found  and  written  to  the  output  catalog  As  with 
body  waves,  a  noise  estimate  is  found  for  the  background  and  a  value  added  tc 
the  amplitude  measurement 

2  3  THE  ANLYST  MODULE 

The  primary  function  of  the  ANLYST  module  is  to  apply  carticular 
network  rules  for  assigning  discrimination  scores  to  the  events  generated  by 
EVENTS  based  upon  the  seismograms  seen  at  the  stations  in  the  network  as 
computed  by  SYNWAV  ANLYST  reads  amplitude,  period,  travel  time  and  noise 
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data  from  SYNWAV  output  files  and  computes  additional  quantities  needed  for 
discrimination  scoring.  These  quantities  are  signal-to-noise  ratios,  the  number  of 
stations  with  good  recordings,  the  body  and  surface  wave  magnitudes.  pP  step- 
out,  and  travel  time  residuals  found  from  locating  the  event.  The  travel  time 
residuals  result  from  a  nonlinear  least  squares  inversion  of  the  travel-time  data 
(read  from  tables).  The  location  program  uses  the  actual  source  coordinates 
written  by  SYNWAV  as  starting  values  for  the  iterative  procedure  so  convergence 
is  rapid.  As  required  by  particular  network  procedures  for  assigning  the  uave1 
time  score,  residuals  are  computed  with  the  depth  fixed  at  the  free  surface  anc 
with  the  depth  a  free  parameter  in  the  inversion  Where  applicable  locations  are 
done  including  the  S  phase.  The  program  then  applies  the  parLcu'ar  net/.or» 
rules  and  writes  a  bulletin  with  the  discrimination  scores 

2  4  THE  DISPLA  MODULE 

The  final  module,  which  was  demonstrated  at  the  Ce'"ter  *zr  S--  .  •• 
Studies  m  November  1983.  allows  the  user  to  1st  arc  pet  ”~e  v-;..  ‘s  .  •  •  • 
analysis  described  above  or  the  obse-vations  of  an  cperat  •  __  -•/.  ■- 
program  has  a  sophisticated  command  processcr  •  . 

specification  of  the  type  of  disp'ays  req^.  red  I"e  , 

divided  mto  two  groups 

1  evert  detection  arc  d  s;n--  •  •  «.  •• 

2  seismic  ty  caramels 

A  th.rd  type  c‘  plot  at  — a ,  :  :•  •• 

depths  of  md'vidua  events  e 

3  M  versus  m  :  vers..s  ..  -  c  .  *•  .  V  •  •  • 

that  sat  s‘y  a  g  ven  co-ct 
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The  data  on  a  given  plot  are  always  based  on  events  that  meet  one  or 
more  conditions  specified  by  the  user.  A  condition  is  a  restriction  on  as  many  as 
eight  properties  of  an  event  (i.e.,  a  condition  has  eight  “subconditions"): 

T  mb, 

2.  Mg, 

3.  depth, 

4.  geographical  area  in  which  the  event  was  located, 

5.  year  in  which  the  event  occurred, 

6.  “identification  status"  assigned  to  the  event, 

7.  discrimination  criterion  on  which  the  status  is  based, 

8.  type  of  event  (mainshock  versus  aftershock). 
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III.  USERS’  MANUAL  FOR  NICE 

This  section  is  a  users'  manual  for  the  NICE  program.  There  are  actually 
four  programs:  events,  synwav.  anlyst  and  displa  Each  is  documented  in 
UNIX  ‘man’'  format. 

The  programs  are  linked  together  through  files.  There  are  10  files  and 
they  serve  as  the  principal  inputs  and  outputs  for  the  programs  The  manual  in 
this  section  includes  man  pages”  for  these  files,  in  addition  to  the  programs 
Particular  file  names  referenced  in  this  section  are  those  used  in  the  testing  and 
initial  calibration  runs  of  the  program.  They  are  included  as  illustration 

The  methods  and  algorithms  for  the  NICE  programs  are  described  in 
Section  II  and  Appendices  A  and  B. 
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PROGRAM 


events  .  .  .  Generate  a  random  seismic  event  catalog. 


SYNOPSIS 


events 


DESCRIPTION 


-a  <  area-list  >  -y  <  year-spec  >  -n  <  num-of-events  > 
-s  <  parameter-file  >  »o  <  output-file  > 


Program  events  generates  a  pseudo-random  catalog  of  earthquakes, 
distributed  in  space,  time  and  magnitude  according  to  a  given  set  of  input 
seismicity  parameters.  The  command-line  arguments  are  defined  below  (Each 
of  these  arguments  is  optional,  and  takes  a  default  value  when  omitted.) 

<  area-list  > . 

The  argument  <  area-list  >  is  a  list  of  geographic  areas  for  which  events 
are  to  be  generated.  Each  entry  in  the  list  refers  to  one  of  five  predefined  areas: 

hkA  .  .  .  Hindu  Kush  (Subarea  A) 
hkB  .  .  .  Hindu  Kush  (Subarea  B) 
hkC  .  .  .  Hindu  Kush  (Subarea  C) 
hkD  .  .  .  Hindu  Kush  (Subarea  D) 
hk  ...  Hindu  Kush  (All  subareas) 

A  given  area  is  referred  to  by  the  short  name  shown  above  (e  g  ,  hkB)  or 
by  a  number  from  1  to  5.  which  selects  according  to  the  above  order  (e  g  .  hkB 
is  2).  Entries  in  --  area-list  >  must  be  separated  with  whitespace:  e  g  .  -a  i  2  3 

The  default  value  of  <  area-list  >  is  T 

<  year-spec  > . 

•  year-spec  ■  is  a  specification  of  which  years  events  are  to  be 
generated  for  There  are  two  acceptable  forms: 


m 


v  v  O  O  v’  I.*  c*  •  •  »  •  .  w  *  »  ‘  .*  . 
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(1)  Two  numbers  separated  by  a  colon  (no 
whitespace),  which  indicate  the  first  and  last 
years  desired;  e.g.,  ‘-y  73:77'. 

(2)  A  single  number,  giving  the  total  number  of  years 
desired.  In  this  case,  the  years  are  numbered 
from  1 ;  e.g,,  ‘-y  4’,  which  is  equivalent  to  -y  1  A'. 

The  default  value  of  <  year-spec  >  isT. 

<  num-of-events  > . 

This  argument  is  the  total  number  of  events  to  be  generated.  The 
distribution  of  seismicity  with  respect  to  time  ts  adjusted  such  that  <  num-of- 
events  >  are  generated  over  the  time  period  indicated  with  <  year-spec  >  . 
Event  origin  times  are  incremented  by  the  recurrence  intervals  At,  given  in 
Appendix  A. 3. 

The  default  value  of  <  num-of-events  >  is  0. 

<  parameter-file  > . 

This  is  a  filename  "prefix”  for  the  input  seismicity  parameter  files.  There 
must  be  one  such  file  for  each  geographic  area  specified  by  <  area-list  >  The 
program  expects  each  such  file  to  be  named  as  follows: 

<  pardir  >  /  <  parameter-file  >  .  <  area-name  > 

<  pardir>  is  the  name  of  a  directory  which  the  user  must  define  with  a  Unix 
environment  variable  called  PARDIR.  This  may  be  accomplished  with  the  Unix 
command  setenv  as  follows: 

setenv  PARDIR  <  pardir  > 

<  parameter-file  ■»  is  supplied  as  an  argument  of  the  -s  flag.  ■  area-name  •  is 
the  name  of  the  area  selected  through  the  -a  flag 
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For  example,  if  you  use 

setenv  PARDIR  /usr/ljacks/bunyon 
events  -s  myparms  -a  hkA  hkC 

then  events  will  expect  seismicity  parameters  to  be  in  the  two  files: 

/usr/ljacks/bunyon/myparms.hkA 

/usr/ljacks/bunyon/myparms.hkC 

The  content  and  format  of  the  seismicity  parameter  files  is  described  in 
the  manual  entry  “Seismicity  File  “ 

The  default  value  of  <  parameter-file  >  is  ‘par’. 

<  output-file  > . 

Filename  for  the  output  event  catalog.  The  full  file  pathname  is 
<  catdir  >  /  <  output-file  > 

where  <  catdir  >  is  the  name  of  a  directory  which  the  user  must  define  with  a 
Unix  environment  variable  called  CATDIR.  this  may  be  accomplished  with  the 
Unix  command  setenv  as  follows: 

setenv  CATDIR  <  catdir  > 

For  example,  if  you  use 

setenv  PARDIR  /usr/ljacks/bunyon 
events  -o  mycatalog 

then  events  will  write  the  event  catalog  it  generates  to  the  file 

usr/ljacks  bunyon  mycatalog 
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If  the  -o  flag  is  omitted,  the  event  catalog  is  written  to  '  standard  output 
i.e.,  the  CRT  screen  or  a  redirected  output  file. 

The  content  and  format  of  the  output  event  catalog  is  described  unoer 
manual  entry  "Event  Catalog." 


PROGRAM 


synwav  .  ...  Generate  a  synthetic  arrival  bulletin. 

SYNOPSIS 

synwav 

DESCRIPTION 

Program  synwav  generates  a  synthetic  arrival  bulletin  for  a  given 
seismic  network  and  set  of  earthquakes.  For  each  earthquake,  the  bulletin 
contains  a  report  of  the  various  seismic  wave  arrivals  (Rayleigh  P  pp  or  S) 
detected  by  each  station  of  the  network.  The  report  includes  amplitude  arrival 
time  and  signal-to-noise  information  about  each  arrival 

The  arrival  bulletin  is  output  to  the  following  file 

/USERS  geop/barker  nice  files  stnbul/hkC 

See  manual  entry  "Arrival  Bulletin"  for  a  description  of  the  file  content  and 
format 

All  input  to  synwav  is  from  files  The  type  and  name  of  each  input  file 

follows 


Station  File 

This  file  contains  the  names  and  locations  of  the  network  stations  (See 
manual  entry  Station  File  ")  The  filename  is 

USERS  geop  barker  nice  files  stnfiles  john 


Event  Cataioc 


This  file  contains  the  catalog  of  events  for  wnicn  a'rva  s  are  tc  be 
generated  It  includes  the  location,  depth  s^ze  and  s  i c  parameters  cf  eac*"’ 
event  (See  manual  entry  Event  Catalog  The  file  name  <s 


USERS  geop  barker  nice  files  events  hkC 
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Earth  Structure. 

This  file  contains  the  earth  structure  used  for  calculating  surface-wave 
and  body-wave  responses.  (See  manual  entry  "Structure  File"  for  more 
information.)  The  filename  is 

/USERS/geop/barker/nice/filesstr/ekz 

Path  File 

This  file  contains  surface-wave  response  information  for  the  assumed 
earth  structure  model.  (See  manual  entry  "Path  File"  for  further  information.) 
The  filename  is 

/USERS  geop  barker  nice,  files  eigfiles  ekz3 

Body  Wave  Table  Files. 

There  is  one  body  wave  table  file  for  eacn  body  wave  to  be  modeled 
(e  g  .  P)  The  file  contains  a  tabulation  of  travel  time  and  amplitude  as  a  function 
of  epicentral  distance  and  event  focal  depth.  (See  manual  entry  Table  File.'') 
The  filename  containing  the  body  wave  table  for  the  wave  type  •  wave  is  given 


/USERS  geop  barker  nice  files  bwtables  tab  wave 
*•  wave  •  is  one  of  P,  pP  or  S 
Body  Wave  Reference  Files 

These  files  supplement  the  body  wave  table  files  There  is  one  reference 
file  for  each  type  of  body  wave  It  contains  information  regarding  the  reference 
values  of  moment,  position  on  the  radiation  pattern  and  reciever  amphhcat  on 
used  m  generating  the  amplitude  table  contained  m  the  correspond  ng  body 
wave  table  file  (For  further  information,  see  the  manual  entry  Re'erence  F  ie  j 
The  filename  containing  the  body  wave  reference  data  ‘or  me  wave  type 
•  wave  ■  is  given  as 

USERS  geop  barker  nice  files  reftab  ref.  wave 
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Long  Period  Instrument  File 

This  file  contains  the  frequency  response  of  the  long-penoo  instrument 
at  each  station.  It  is  assumed  to  be  station-independent.  (See  manual  entry 
"Instrument  File.")  The  filename  for  the  long-period  response  is 

/USERS/geop/barker/nice/files'inst  kono.lpz 

Short  Period  Instrument  File 

This  file  contains  the  frequency  response  of  the  short-period  instrument 
at  each  station.  It  is  assumed  to  be  station-dependent  (See  manual  entry 
"Instrument  File  ")  The  filename  for  the  short-period  response  is 

/USERS/geop/barker/nice  files  inst  kono.spz 
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PROGRAM 

anlyst  .  Generate  a  synthetic  discrimination  bulletin 


SYNOPSIS 

anlyst 

DESCRIPTION 

Program  anlyst  generates  a  synthetic  discrimination  bulletin  for  a  given 
seismic  network  and  set  of  earthquakes.  The  bulletin  is  obtained  by  applying 
detection,  location  and  discrimination  algorithms  to  the  data  in  an  arrival  bulletin. 
(See  manual  entry  'Arrival  Bulletin".)  For  each  event  reported  in  the  arrival 
bulletin,  anlyst  produces  a  record  for  the  discrimination  bulletin  which  contains 
an  estimate  of  the  event's  origin  time,  location  and  magnitudes  as  well  as  a  set 
of  "scores"  which  summarize  the  results  of  discrimination  analysis. 

The  discrimination  bulletin  produced  by  anlyst  is  output  to  the  following 
file: 


USERSgeop/barker/nice/files/aftac/hkC 

See  manual  entry  "Discrimination  Bulletin"  for  a  description  of  the  file  content 
and  format 

All  input  to  anlyst  is  from  files.  The  type  and  name  of  each  input  file 

follows. 

Station  File. 

This  file  contains  the  names  and  locations  of  the  network  stations  (See 
manual  entry  Station  File".)  The  filename  is 

USERS  geop  barker  nice  files  stnfiles  john 
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Arrival  Bulletin 


This  file  contains  the  reports  of  seismic  wave  arrivals  (e  g  .  P.  oP.  S 
Rayleigh)  at  the  network  stations  from  various  events  (See  manual  entry 
"Arrival  Bulletin".)  The  filename  is: 

/USERS/geop/barker/nice/files  stnbul  hkC 

Body  Wave  Table  Files. 

anlyst  requires  one  body  wave  table  file  for  each  of  the  body  waves:  P. 
pP.  S.  Each  file  contains  a  tabulation  of  travel  time  and  amplitude  as  a  function 
of  epicentral  distance  and  event  focal  depth.  (See  manual  entry  Table  File  '.) 
The  filename  containing  the  body  wave  table  for  the  wave  type  •  wave  ■  is 
assumed  to  be 

/USERS/geop/barkernice/files/bwtables'tab  <  wave  * 
where  <  wave  >  is  one  of  P,  pP  or  S 
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PROGRAM 

displa  .  .  Query  and  display  a  discrimination  bulletin 


SYNOPSIS 

displa 


DESCRIPTION 

displa  is  an  interactive  program  for  displaying  the  oata  in  a 
discrimination  bulletin  The  bulletin  contains  a  summary  of  discrimination  results 
for  a  given  set  of  events.  Displays  take  the  form  of  plots  or  tables  showing 
various  breakdowns  of  these  results  (Plots  require  a  Tektronix  4014  terminal  ) 

The  content  of  a  given  display  is  governed  by  two  types  of  information 
specified  by  the  user.  The  first  type  of  information  is  a  set  of  conditions 
determining  a  subset  of  the  entire  set  of  events  m  the  bulletin  For  example,  the 
user  might  specify  the  subset  of  events  which  satisfy  the  two  conditions  (i)  the 
event  has  been  correctly  identified  by  the  network  and  (2)  its  focal  depth  is  less 
than  50  km  More  than  one  such  subset  might  be  involved  m  a  given  display 

The  second  type  of  information  specifies  the  event  properties  to  be 
displayed  for  those  events  in  the  specified  suosetis)  Examples  cf  such 
properties  are  the  distributions  of  the  events  witn  resDect  to  magnituce  and 
depth 

Now  we  define  these  two  types  of  .nformation  m  oeta 
Event  Subsets 

A  given  subset  of  events  'S  oefmea  by  conditions  cn  tne  foiicwmp  e  arn 
event  properties 

Ci  The  fype  of  the  event  main  shocx  ahe-sioo  0'  a1'  event 
which  nas  been  ignored  for  analysis 


C2  The  class  of  the  event,  as  determined  by  its  discrimination 
scoring  results. 

C3.  The  discriminants  whose  scores  are  "credited''  with 
determining  the  class  of  the  event. 

C4.  The  geographical  area  where  the  event  has  been  locatec. 
The  possible  areas  have  been  predefined.  Arbitrary 
conditions  on  latitude  and  longitude  are  not  allowed. 

C5  The  year  in  which  the  event  occurred. 

C6.  The  body  wave  magnitude,  determined  for  the  event 

C7.  The  surface  wave  magnitude.  M  determined  for  the 
event. 

C8  The  estimated  depth,  z.  of  the  event. 


A  condition  on  the  latter  three  quantities  (m^,  Mg.  z)  consists  of  a  minimum  and 
maximum  value. 

The  actual  commands  for  specifying  event  subsets  through  the  above 
conditions  are  described  later. 

Functions  to  Display 

One  may  display  the  following  event  properties  mfa.  z  anc  event 
counts  and  proportions  as  a  function  of  these  three  variables  These  cuantTes 
can  be  displayed  in  many  different  ways 

To  describe  the  possible  displays,  let  x  and  y  each  denote  o~e  of  me 
variables  m^.  Mg  or  z.  Let  A  and  B  denote  two  subsets  of  events  whicr  tre  user 
has  defined  through  a  list  of  conditions  The  data  from  the  events  -n  A  ana  B 
may  be  displayed  as  the  following  functions 

p 

Ft  N^(x).  the  cumulative  count  of  events  m  A  with  respect  tc 

x  i  e  .  N’tfXp.)  is  the  number  of  events  in  A  with  x  >  x« 

A'  0  0 


I 


F2.  NA(x):  the  incremental  count  of  events  in  A  with  respect  to 
x  Given  a  specified  increment  Ax.  NL(Xq)  is  the  number  of 
events  such  that  xq  <  x  <  Xq  -  Ax.  Thus, 

Na(x)  =  N$=  (x)  -  NCa(x  -  Ax) 

F3.  PA)  g(x):  the  cumulative  ratio  of  the  number  of  events  m 
A*B  to  the  number  in  B.  A*B  denotes  the  intersection  of 
subsets  A  and  B.  Note  that 

PAlB(x)  =  nCA’BW  N°B(x) 

Observe  that  N  .«□  is  the  number  o'  B  events  wn  cn  are 

Ad  r 

also  A  events.  Thus  A’B  is  a  subset  of  B  anc  P  A  (  □  must 
be  between  0  and  1 

F4.  PA  g(x):  the  incremental  ratio  of  the  number  of  events  n 
A*B  to  the  number  m  B 

PAlB(x)  =  NA’B(x)  NBfx) 

F5  y(x)  distribution  of  events  (from  subset  A)  m  the  x-y  p'ane. 
e  g  .  Ms(mb)  defines  the  classical  versus  m^  plot 

Commands  in  displa 

After  executing  displa  the  user  is  quer  ed  for  the  name  o'  a  file 
containing  a  discrimination  bulletin  She  will  then  be  asKec  tc  instate  whether 
this  file  is  formatted  or  unformatted  (If  the  file  has  been  generateo  by  anlyst.  it 
is  formatted .)  After  this,  the  user  operates  displa  through  a  set  of  commands 

There  are  four  types  of  commands  recognized  by  displa 

1  Commands  to  define  conditions  on  evert  subsets 

2  Commands  to  display  data  from  the  discrimmat  or  bulletin 

3  Commands  to  set  parameters  for  displays 


4  Utility  commands 
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Condition  Commands  (A  Conditions). 

Earlier  we  defined  various  quantities  one  could  disolay  for  a  subset  of 
events  we  called  A  or  for  two  subsets  A  and  B;  e.g.,  NA,  PA  (  g  At  any  given 
time  during  a  run  of  displa.  A  and  B  are  defined  in  terms  of  conditions  on  eight 
variables  (see  C1-C8  above).  In  this  section  we  discuss  commands  for  the  user 
to  specify  “A  conditions";  i.e.,  conditions  which  determine  the  event  subset  A. 

When  the  user  initiates  a  run.  a  set  of  default  conditions  are  invoked. 
These  define  the  subset  A  of  events  to  be  all  inclusive:  i.e..  A  is  the  entire  set  of 
events  in  the  bulletin.  To  redefine  A,  the  user  may  specify  her  own  conditions 
with  one  or  more  commands  of  the  form: 


<cond-var>  =  <  A-cond-value  > 

(The  equals  sign  (  =  )  is  optional.)  For  <  cond-var  >  one  uses  one  of  the 
following  keywords: 


ty 

cl 

di 

ar 

yr 

mbr 

Msr 

zr 


to  specify  event  types 
to  specify  event  classes 

to  specify  discrimination  criteria  for  event  classes 

(explained  below) 

to  specify  areas 

to  specify  years 

to  specify  m b  ranges 

to  specify  M  ranges 

o 

to  specify  z  ranges 


At  any  given  time,  a  condition  command  for  a  given  variable 
automatically  overrides  any  previous  definition  for  that  variable,  whether  it  was  a 
default  or  user  specified.  Conditions  on  other  variables,  however,  are  not 
affected 

There  are  two  styles  for  <  A-cond-value  • .  depending  on  which  variable 
one  is  specifying  For  mbr,  Msr  and  zr.  *-  cond-var  •  specifies  a  minimum  and 
maximum  value: 
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<  cond-var  •  =  <  mm-vaiue  --  -max-value  - 

This  command  restricts  the  subset  A  to  include  only  events  navmg 

<-  mm-value  >  <  cond-var  >  -  •-  max-value  -- 

For  example,  to  include  events  with  Mg  between  5  0  and  6  5  (not  including  6  5). 
one  should  use  the  command 

Msr  -  5:6  5 

The  second  style  for  ■  A-cond-vaiue  •  selects  one  or  more  items  from  a 
fixed  menu  of  choices  This  style  applies  to  ty  cl  di  ar  and  yr  The  simplest 
form  is  to  specify  a  single  item  i  e 

<  cond-var  •  =  -  item 
For  example. 

yr  =  74 

To  select  two  or  more  items  one  uses 

•  cond-var  •  -  •  item1  •  -  •  item  -  -  -  -  tem^ 

In  this  way.  an  event  is  in  subset  A  if  it  matches  any  of  toe  items  g1.en  Fc' 
example,  to  select  events  m  the  president^  eiectioo  years  form  me  s 

yr  -  72  -  76 

In  general,  one  may  use  any  iogica  expression  mvc  v  ng  :-ie  re-u  'ter"s 
and  the  logical  operators  -  (or)  *  (and;  and  (not  Tnus  72  76  '"er-s 
"events  in  the  years  72  or  76"  As  another  example  to  se  ect  e.e-ts  '■  m  a 
years  except  1973  (i  e  .  not  73  )  one  would  use 
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yr  =  73 

Let  us  now  define  the  complete  menu  for  each  condition  variable 

Mi.  Types  ma  af  ig 

ma  The  event  is  a  main  shock. 

af  .  .  The  event  is  an  after  shock. 

ig  .  The  has  been  ignored  (not  analyzed). 

M2  Classes  Cl  C2  C3  C4  c5  C6  c7  c8 

The  classes  are  defined  in  terms  of  numerical 
discrimination  scores  as  follows.  The  following  definitions 
are  hierarchical  i  e  .  the  first  test  which  the  scores  satisfy 
determines  the  class. 


cl 

one  or  more  1  s 

c2 

two  or  more  2's. 

c3 

one  2 

c6 

two  or  more  4's  or 
one  4  and  one  or  more  3's 

c5 

one  4.  or 

two  or  more  3's 

c4 

one  3 

c7 

one  or  more  5  s:  or 
one  or  more  6's 

c8 

no  scores  given 

M3  Discriminants  dl  d2  d3  d4  d5  d6 

M4  Areas  kkn  kks  es  hk  cau 

These  are  abbreviations  for  geographic  areas,  defined  as 
follows 
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cau 


Kurile-Kamchatka  north 

(longitude  >  140*E.  latitude  >  48*N) 

Kurile-Kamchatka  south 

(longitude  >  140*E,  latitude  <  48*N) 

Eastern  Siberia  (85*  <  longitude  < 

1 40  *  E) 

Hindu  Kush  (85*  E  <  longitude  <  60* E) 
Caucasus  (longitude  <  60*  E) 


M5  Years:  73  74  75  76  77 

A  condition  on  discriminants  is  actually  an  elaboration  on  the 
condition  on  class.  Each  class  requires  that  certain  discrimination  score 
values  occur.  Discriminants  that  have  the  proper  score  values  are 
“credited"  with  the  event  having  been  assigned  to  its  class.  The 
condition  on  discriminants  thus  restricts  the  subset  of  events  to  include 
only  those  event  for  which  the  specified  discriminants  receive  credit  for 
the  event's  class.  For  example,  events  which  are  in  class  1  by  virtue  of 
the  discriminant  d4  are  selected  with 

cl  =  cl 
di  =  d4 

Given  the  definitions  of  the  classes,  this  is  equivalent  to  the  set  of  events 
for  which  discriminant  d4  received  a  score  of  1 .  As  another  example,  to 
select  events  for  which  both  d2  and  d4  receive  credit  for  either  class  1  or 
class  2,  one  uses 

Cl  =  Cl  -  C2 
di  =  d2  *  d4 

This  is  satisfied  by  either  (a)  events  in  class  i  for  which  d2  and  d4  both 
have  a  score  of  i ,  or  (b)  events  in  class  2  for  which  o2  and  d4  both  have 
a  score  of  2 

Note  that  a  simple  conjunction  of  items  using  *  or'y  makes  sense  m  the 
case  of  conditions  on  discriminants  That  is.  two  or  more  discriminants  car 
simultaneously  receive  credit  for  an  event's  class  However,  an  event  cannot,  for 
example,  be  simultaneously  in  two  or  more  areas  or  in  two  or  more  years  Thus 
for  example  ar  =  kkn  *  kks  yields  the  empty  set  of  events 


As  a  convenience  for  experienced  users,  the  menu  items  *or  each 
condition  variable  are  given  a  fixed  order,  and  the  user  may  use  the  mm-max" 
style  to  select  a  set  of  contiguous  items  from  a  menu  For  example,  tne  ordering 
of  years  is  chronological,  so 

yr  =  73  78 

is  equivalent  to 

yr  =  73  -  74  *  75  -  76  -  77  *  78 

As  a  further  convenience,  the  menu  items  are  numbered  and  the  user  may 
select  items  by  number  (as  long  as  she  has  figured  out  the  numbering).  For 
example, 

ar  =  2 

is  equivalent  to 

ar  =  kks 

As  a  final  convenience,  the  alias  all  has  been  defined  to  denote  the  "all- 
inclusive"  condition  for  any  variable  For  example. 

ar  =  all 

denotes  ar  =  kkn  -  kks  -  es  -  hk  -  cau  ,  and 

mbr  =  all 


denotes  'mbr  =  -1000:1000'  (i.e  .  or  some  other  suitably  large  magnitude 
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To  learn  the  current  setting  of  a  given  condition  variable,  tyoe 

<  cond-var  >  =  ? 
eg.,  ‘ar  =  ?’. 

Condition  Commands  (B  Conditions). 

When  the  user  initiates  displa.  the  subset  B  is  undefined  To  define  B 
an  extended  syntax  of  condition  commands  is  recognized: 

<  cond-var  >  =  <  A-cond-value  >  '  <  B-cond-vaiue  > 

<  A-cond-value  >  specifies  a  condition  on  subset  A.  as  described  above  <  B- 
cond-value  >  specifies  a  condition  with  respect  to  the  same  variable  on  event 
subset  B. 

Once  a  B  condition  has  been  explicitly  defined  for  one  variable,  displa 
automatically  invokes  a  default  value  for  the  B  conditions  of  the  other  variables 
The  default  is  the  corresponding  A  condition.  Thus. 

<  cond-var  >  =  <  A-cond-value  > 

is  the  same  as 


<  cond-var  >  =  <  A-cond-value '»  •-  A-cond-vaiue  ■ 

For  example,  the  following  sequence  defines  A  and  B  to  have  events  m  the  sane 
two  areas,  but  with  different  classes 

ar  =  kkn  *  kks 
cl  =  cl  /  c1:c7 

In  this  example.  P^(  g  would  be  the  fraction  of  events  ir  the  Kunie-KamchatKa 
areas  which  are  m  class  1.  among  those  from  the  same  two  areas  which  are  m 
classes  i  through  7 
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Condition  Commands  (Multiple  Conditions, 

It  is.possibie  to  define  multiple  versiors  of  tne  suoset  A  Arne-  m  s  s 
done,  displa  automatically  generates  multiple  functions  on  a  smgie  d  so  ay  cms 
facilitating  the  comparison  of  data  For  example  suppose  one  warned  tc 
compare  data  from  events  in  three  different  geographic  areas  One  would  men 
define  three  subsets  A1 ,  A^  and  A^  which  satisfied  a  common  set  of  conations 
except  for  the  condition  on  area  A  display  o*  a  given  function  would  tr,er 
automatically  show  three  versions  one  for  each  area  If  one  reouesteo  a  disc  ay 
of  for  example  the  three  functions  N A ., ( m p )  NAO(mb)  anc  m^ 

would  appear  on  the  same  display 

Multiple  subsets  A  may  differ  with  respect  tc  only  one  concihon  var  ac  e 
The  command  syntax  for  defining  multiple  versions  of  tne  desired  variaDie  is 

-  cond-var -•  =  •  A1 -cond-value  •.  ■  A^-cond-vaiue  ■ 

Thus,  multiple  condition  values  are  separated  by  commas  Each  value  follows 
the  rules  outlined  earlier  As  an  example,  suppose  one  wants  to  make  an  area- 
by-area  comparison  of  class  1  events  from  19?3  Then  one  wou'd  de'^e 
multiple  A  subsets  with 

ar  =  kkn  kks  es  hk  cau 

cl  =  1 

yr  =  73 

These  commands  imply  that  A(.  i  =  1 .  .5  are  definec  by  the  same  ccmd  t  ons 

on  class  and  year,  but  different  conditions  on  area 

The  reader  should  be  careful  to  distinguish  between  come  ~  rg  .te-^s 
with  '  -  '  and  defining  multiple  items  with  Tne  -  syntax  merges  two  te~s 
mto  one  while  '  separates  individual  cases  Tc  emcnas  ze  tne  di‘ferer'ce  rote 
the  following  example  wnich  uses  both  devices  Suppose  we  wis^  tc  scream 
events  from  two  geographic  regions  tne  Kun  e-Kamcnatka  'egio*'  ve'Sus  me 
rest  of  the  Eurasia  We  would  then  cef  ne  two  suosets  of  events  cy  sees  ‘v  '~c. 
multiple  conditions  on  area 
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ar  =  kkn  -  kks  es  -  hk  -  cau 

Thus,  the  first  region  (subset  A1 )  combines  kkn  and  kks.  The  second  subset 
combines  events  from  three  areas 

Multiple  B  subsets  may  also  be  defined,  but  only  implicitly:  i  e  ,  multiple 
items  may  not  appear  behind  a  Multiple  versions  of  B  are  implicitly  defined 
for  a  variable  when  multiple  A  conditions  have  been  given  for  the  variable  but  no 
B  condition  has  been  given  For  example,  suppose  we  set  every  conation 
variable  to  all.  except  for  the  following 

ar  =  kkn  kks  es 
cl  =  cl  all 

This  defines  three  subsets  A1 ,  A2  A^  and  three  subsets  B1 .  B2  The  A  s 
contain  events  m  class  1.  while  the  B  s  contain  events  in  any  class  (i  e  .  from  all 
eight  classes)  A1  and  B1  contain  events  m  area  kkn.  A2  and  B2  for  kks.  and  so 
on  Given  these  subsets,  if  one  were  to  ask  for  a  display  of  N^m^).  she  would 
get  three  functions  N^m^).  i  =  i  ,2.2  Each  is  the  number  of  class  1  events  in 
its  respective  area  If  one  were  to  ask  to  display  PA|  then  she  would  get 

three  functions  PAi  i  B/mb^'  1  =  1,2  3  PAi  I  Bi  (mb^  lS  the  fractl°n  °f  class  1 
events  m  kkn  among  all  events  in  kkn  And  so  on 

Now  compare  this  to  the  following 

ar  kkn  kks  es  all 
cl  -  cl 

This  defines  three  subsets  A1  A.,  A^  and  one  subset  B  The  A  s  are  as  before 
the  number  of  class  i  events  m  each  of  three  areas  B  is  number  of  class  i 
events  m  ail  areas  combined  (To  understand  this  recall  that  cl  cl  is  the 
same  as  cl  cl  cl  )  Now  N^im^  oehnes  the  same  three  ‘unctions  as 
cefme  PAlB^b)  though  refers  to  PA1 , B,mb;  PA2,Btmb:  ana  PA3|BitrtD) 
PA1 1  g  is  the  fraction  of  class  i  events  m  area  kkn  among  class  i  events  m  all 
areas  combined 
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As  an  additional  convenience  for  defining  multiple  conditions  tnere  s  an 
alias  1  which  selects  each  and  every  menu  value  for  a  conoition  variable  For 
example. 

ar  =  ! 

is  equivalent  to  ar  =  kkn  kks,  es,  hk.  cau. 

Display  Commands 

There  are  two  commands  for  displaying  data 

plot  -  function  > 
list  <-  function  • 

plot  produces  a  standard  x-y  plot  containing  one  or  more  curves  list  produces 
a  table  with  one  column  for  x-axis  data,  followed  by  one  or  more  columns  for  y- 
axis  data.  The  data  whicn  are  displayed  are  defined  by  the  argument 
<  function  -• . 

The  argument  <  function  >  takes  the  general  form  <  y  >  ( •  x  ' ),  where 
•:  y  >  denotes  a  dependent  variable  and  •-  x  •  denotes  an  independent  variable 
The  choices  for  x  >  are: 


m.  M.  z 
b  s 


with  obvious  meaning  The  choices  for  <y>  are 


nc 

n 

pc 

P 


cumulative  count  of  events  in  subset  A.  Nc^ 
incremental  count  of  events  in  subset  A. 
cumulative  ratio  of  events  in  subsets  A  and  B.  Pc^ ,  g 
incremental  ratio  of  events  in  subsets  A  and  B.  PA  j  g 


For  plots  (but  not  lists)  *  y  ■  may  also  be 


m.  z 
b  s 
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When  A  and ''or  B  are  multiply  defined,  multiple  plots  or  tabulations  are 
automatically  produced  for  A1 .  .  and  B  or  B1  B2  as  describee  -n  tne 

previous  section. 

Examples  of  display  commands  are 

plot  nc  (mb) 
list  p  (Ms) 
plot  Ms  (mb) 

As  a  convenience,  displa  allows  the  combination  of  display  and 
condition  commands,  with  the  forms 

plot  <  function  >  <  condition1  >  <  condition2  • 
list  <  function  >  <  condition^  >  <  condition^  •  ■ 

Each  of  the  optional  arguments  <  conditio^  >  ....  is  a  condition  command  as 
defined  above  It  behaves  exactly  as  if  it  had  been  issued  prior  to  the  display 
command,  including  the  fact  that  it  alters  the  definitions  of  A  and  B  for  future 
commands  (until  overridden)  Thus 

plot  pc  (mb)  ar  =  1  cl  =  cl  -  all 

is  equivalent  to 


cl  =  cl  /  all 
ar  -  i 
plot  pc  (mb) 

Parameter  Commands. 

Parameter  commands  set  certain  parameters  for  plots  and  lists  There 
are  two  types  of  such  parameters 
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1.  Plot  table  limits'  i.e  ,  the  minimum  ana  maximum  values  of 
the  x-axis  variables  for  which  data  are  displayea 

2.  Plot/table  increments:  i.e..  the  sampling  increment  for  the  x- 
axis  variable,  determining  the  spacing  between  displayed 
data. 

There  are  three  possible  x-axis  variables  (m^.  Mg,  z).  so  there  are  six  commands 
for  setting  display  parameters: 


mblim 

Msiim 

zlim 

and 

mbinc 

Msinc 

zinc 


=  <  mm-value  > 
=  <  mm-value  > 
=  <  mm-value  > 


=  <  value  > 
=  <  value  > 
=  <  value  > 


<  max-value  > 

<  max-value  > 

<  max-value  > 


To  learn  the  current  setting  of  these  parameters,  type: 
<  parameter  >  =  ? 


eg,  mbinc  -  ?’. 

Utility  Commands 

The  following  additional  commands  are  recognized  by  displa 

help  print  out  some  information  which  might  help  the 

user  select  his  next  command 


quit 


terminate  this  session  with  displa 


FILE  TYPE 


Arrival  Bulletin  report  of  seismic  arrival  mformat  on  ,rcn 

a  set  of  events  at  a  set  of  stations 


CREATED  BY 

Program  synwav 


USED  BY 

Program  anlyst 


DESCRIPTION 

The  arrival  bulletin  contains  the  arrivals  reported  by  a  network  for  each 
of  several  seismic  events.  The  arrival  information  includes  arrival  times, 
amplitudes,  periods,  and  noise  levels  for  various  seismic  phases 

The  arrival  bulletin  is  an  ASCII  file  consisting  of  records  separated  by 
new-lmes  The  records  are  in  groups,  each  group  corresponding  to  one  event 
and  its  arrivals  Within  a  group  is  a  set  of  records  for  each  station,  each  set 
containing  information  about  the  Rayleigh  wave  arrival  and  various  body  wave 
arrivals 

More  specifically,  tne  bulletin  is  organized  as  follows 

(For  each  event 

Event  header  record 
(For  each  station 
Station  record 
Surface-wave  record 
Body-waves  header  record 
(For  each  type  of  body  wave 
Body -wave  record  i 
Body-wave  record  2 


The  formats  of  the  various  records  follow 


Event  Header  Record 


Fortran 

Columns  Format  Data  Description 

i-5i  v  Event  information  (time,  location,  depth 

moment  focal  mechanism) 

52-55  i4  Number  of  stations  recording  event 

(a^  See  manual  page  on  Event  Catalog " 

Station  Record 

Fortran 

Columns  Format  Data  Description 


1-10 

a 

Station  identifier 

11-20 

i pel 0  2 

Event-station  distance  (degrees) 

21-30 

ipeiO 2 

Event-station  azimuth  ''decrees  c  v 
north; 

31-40 

i pel 0  2 

Event-station  caci* -azimuth  cegre' 
from  ncrth.j 

Body-Wave  Record  (First  of  Two). 
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Fortran 

Columns  Format 


Data  Description 


1-10 

a 

Wave  type  (P\  p 

11-40 

Onset  arrival  time 

11-14 

i4 

Year 

15-18 

i4 

Month 

19-22 

i4 

Day 

23-26 

i4 

Hour 

27-30 

i4 

Minute 

31-40 

f  1 0.2 

Seconds 

Body-Wave  Record  (Second  of  Two). 


Fortran 

Columns  Format 


Data  Description 
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reportcf  o  c  £  t*  o  r  a°c 
discrimmatior  results  fcr  s  set 
of  earthquakes 

CREATED  BY 

Program  anlyst 

USED  BY 

Program  displa 


FILE  TYPE 

Discrimination  Bulletin 


DESCRIPTION 

The  discrimination  bulletin  is  an  ASCII  file  consisting  of  records 
separated  by  new-lines.  Each  record  corresponds  to  one  earthquake  The 
format  of  a  record  is  as  follows: 


Fortran 

Columns  Format  Data  Description 


1-5 

i5 

Event  number 

6-8 

i3 

Year  number  of  origin  time  (earliest  year 
m  bulletin  is  1 ,  secona  earliest  is  2  etc  ) 

10-12 

i3 

Area  number  of  location 

14-16 

i3 

Focal  depth  (kilometers) 

18-20 

i3 

Body  wave  magnitude,  m. 

(hundredths  of  a  magnitude  unit) 

22-24 

i3 

Surface  wave  magnitude  M 
(hundredths  of  a  magnitude^jmt) 

Azimuth  of  major  axis  of  location  error 
ellipse  (deqrees  c  w  from  nodnj 
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30-32 

i3 

Length  of  semi-major  axis  of  location 
error  ell'pse  (kilometers) 

34-36 

i3 

Length  of  semi-minor  axis  of  location 
error  ellipse  (kilometers) 

38-40 

i3 

Event  type  flag 

42-44 

i3 

Score  for  discriminant  1 

46-48 

i3 

Score  for  discriminant  2 

50-52 

i3 

Score  for  discriminant  3 

54-56 

i3 

Score  for  discriminant  4 

58-60 

i3 

Score  for  discriminant  5 

62-64 

i3 

Score  for  discriminant  6 

3-30 
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Event  Catalog  .  .  .  location  and  source  parameters  for  a  set  o* 

earthquakes. 


CREATED  BY 

Program  events 

USED  BY 

Program  synwav 

DESCRIPTION 

The  event  catalog  is  an  ASCII  file  consisting  of  records  separated  Dy 
new-lines.  Each  record  corresponds  to  one  earthquake  The  format  of  a  recora 
is  as  follows 


Fortran 

Columns  Format  Data  Description 


1-6  16 

7-19 

7-8  i2 

9-10  i2 

11-12  i2 

13-14  i2 

15-16  i2 

17-19  i3 

20-32 

20-24  i5 

25  a 

26-31  16 

32  a 


Event  number 

Event  origin  time 

Year 

Month 

Day 

Hour 

Minute 

Seconds 

Event  location 

Latitude  (thousandths  of  degrees  from 
equator) 

Hemisphere  (  n'  or  s') 

Longitude  (thousandths  of  degrees  from 
Greenwich) 

Hemisphere  (  e'  or  w') 


m 


>3 


3-31 


Event  depth  (kilometers) 

Logarithm  of  event  moment  (newton- 
meters) 

Focal  mechanism 

Strike  (tenths  of  degrees  c.w.  from  north) 
Dip  (tenths  of  degrees  below  horizontal) 
Rake  (tenths  of  degrees  c.c  w  from 
horizontal) 
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FILE  TYPE 


Instrument  File  ...  a  seismic  instrument  response 


CREATED  BY 
User. 


USED  BY 


Program  synwav 


DESCRIPTION 

This  file  contains  a  seismic  instrument  response  in  the  form  of  poles  and 
zeros  in  the  complex  frequency  plane. 

Let  s  be  the  complex  Laplace  transform  parameter  (iw).  Then  the 
instrument  response  as  a  function  of  s  is  given  by: 


(s)  =  A  (s  -  z.,)  .  (s  -  zn)  /  [(s  -  p^  (s  -  pm)] 


where  z1 . zp  are  the  zeros  and  p1 . pm  the  poles  of  the  response 

function 

The  instrument  file  is  an  ASCII  file  consisting  of  records  separated  by 
new-lmes  The  data  are  in  Fortran  77  "free"  format  (i  e  .  'read  *')  Therefore, 
arbitrary  white-space  may  be  used  to  separate  the  items  on  each  record 
(including  extra  new-lmes) 

The  contents  of  the  file,  record  by  record,  are  as  follows 


Record  1 

comment 


Record  2: 

<  A  > 

Record  3: 

<  num-of-poles  > 

Record  4: 

<  real-pole..  >  <  imag-pole1  > 

Record  5 

<  real-pole2  >  * ■  imag-pole2  • 

Record  3  -  m: 

<  real-pole  >  *  imag-poiem  * 

Record  4  -  m: 

<  num-of-zeros 

Record  5  -  m 

<  real-zero1  >  <  imag-zero1  • 

Record  6-  m: 

<  real-zero2  <  imag-zero2 

Record  4  -  m  -  n: 

<  real-zerop  >  -  imag-zerop  ■ 

In  the  above,  subscripts  m  and  n  stand  for  the  variables  •  num-of-poies  ■  ana 
•  num-of-zeros  • .  respectively 

The  variable  •  comment  ■  is  a  character  string  of  length  80  or  less  We 
identify  the  other  variables  with  the  variables  m  the  equation  for  I (s )  above,  as 
follows 
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FILE  TYPE 

Path  File  .  .  .  surface  wave  path  and  source  response  functions 

CREATED  BY 
User. 

USED  BY: 

Program  synwav 

DESCRIPTION 

This  file  contains  a  tabulation  of  Rayleigh  wave  path  responses  as  a 
function  of  frequency,  and  source  excitation  eigenfunctions  as  a  function  of 
source  depth  and  frequency 

This  is  a  binary  file  containing  variable  length  records  The  content  of 
each  record  is  as  follows: 

Record  1 : 

<  comment  • 

Record  2: 

num-of-frequenoes  • 

Record  3 

*■-  freq1  •  -  freq2  •  •  freqn  • 

•- phase-ve^  >  •- phase-vei2  ■  •  pnase-ve!n  - 

•  group-vel1  •  -  group-vel2  ■  •  group-velp  • 

*-  a-sub-r1  •  •  a-sub-r2  •  •  a-sub-rn 

ell  ip  1  >  •  eihp2  ■  •  ell  ipn  ■ 

•  7i  "  ^2  "7n  ' 

Record  4 

•  src-depth  -  ■  src-a.  ■  ■  src-/?.  •  src-p.  ■ 
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Record  5: 
<  eiaf 


Record  6: 


111  > 

<e,gf211  >  ■ 

<e'9,n11 

1 12  > 

<  ©iQf 2  2  >  ■ 

■  <e'9fn12 

1 13  > 

<eigl213>  . 

■  •''e'9,n13 

114  > 

<  ^'9^214 

■  <ei9'nl4 

121  > 

<eigf22i  >  • 

-  '  ei9fn2i 

122 

<  ©19(222  >  ■ 

*  eigfn22 

123 

<  e,9f223 

'  ei9'n23 

124  > 

n  6!g^224 '  ’ 

ei9fn24 

1ml  :> 

*  ei9'2m1  ' 

'  e,9'nm1  :’ 

1m2  '' 

'  ei9f2m2  ' 

‘  e'9'nm2  :" 

1m3 

<e'9l2m3  ' 

‘  e'Sfnm3  ;' 

1m4  > 

‘  e'9f2m4  > 

-1  ei9'nm4  s 

Record  4  -  m: 

'e,9'im1 

':e'9fim2 

"  ei9*im3 
'  ei9'lm4 

End-of-file 


In  the  above,  subscript  n  stands  for  <  num-of-frequencies  >  Subscript  m  is  the 
number  of  source  depths  (which  are  counted  as  they  are  read) 

The  meanings  of  the  above  variables  follow. 

■  comment  • 

A  character  string  of  length  80  or  less 
•  num-of-frequencies  - . 

The  number  of  frequency  samples  at  which  the  response  functions  are 
tabulated  (Fortran  type  integer  .) 


The  i  sampling  frequency  (hertz)  (Fortran  type  real 


BE 

j 

A 


•■A  A  .** 


*  v.v; 


phase-vel  >  <group-velf>  <a-sub-r 


These  are  the  Rayleigh  wave  responses  at  the  frequency  ■  freq|  :■ 
These  response  functions  are  those  defined  by  Harkrider  (1S70)  From  left  to 
right,  they  are: 


1.  phase  velocity  (m/s) 

2  group  velocity  (m/s) 

3  Ar  (m/Nt) 

4  ell  i  pticity 
5.  iplkm'1) 

(All  are  Fortran  type  real' ) 

•  src-depth  .■ 

The  jth  sampimg  source  depth  for  eigenfunctions  (Fortran  type  real'.) 

-  src-g[:-  •  src-/9  >  <  src-^  >  <-■  src-^  >  <  src-X)  >  ■-  src-o|  > 

The  earth  properties  at  the  source  depth  <  src-deptlr  •  From  left  to 
right,  they  are: 

1 .  P  wave  velocity  (m/s) 

2.  S  wave  velocity  (m/s) 

3  Density  (kg/m3) 

4.  Shear  modulus  (Pa) 

5  Lame'  parameter  \(Pa) 

6  Poisson's  ratio 

(All  are  Fortran  type  rear  ) 


The  surface-wave  eigenfunction  of  the  kth  type  (k  -  1.2  3.4)  at  the  ilh 
frequency  and  jth  source  depth  (See  Appendix  B  2.)  (Fortran  type  rea'' ) 
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FILE  TYPE 

Reference  File 


CREATED  BY 
User. 

USED  BY 

Program  synwav 

DESCRIPTION 

This  file  is  associated  with  a  particular  body  wave  (e  g  .  P)  and  s  an 
adiunct  to  the  body  wave  table  file  for  the  same  wave  (see  manual  ent'y  Tab’e 
File") 

This  is  an  ASCII  file  It  contains  three  numbers  in  Fortran  77  ,ree  fcrmat 
(i  e  .  read  *').  The  numbers  may  be  separated  by  arbitrary  white-space  b  sens, 
tabs  or  new-lmes)  The  three  numbers  are  the  reference  va  ^es  *zr  wrc"  tre 
amplitude  tables  were  computed  They  may  be  used  to  scaie  tt'e  a^c  t..oe 
tables  The  three  numbers  are  the 

1  The  reference  moment  (N-m) 

2  The  reference  value  on  the  radiation  pattern 

3  The  reference  receiver  amplification  function 


reference  values  for  amphtuce  tables  *cr  a 
padicuiar  body  wave 
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FILE  TYPE 

Seismicity  File  .  .  .  parameters  of  the  seism, city  distribution  n 

a  given  geographic  area 

CREATED  BY 

User. 

USED  BY 

Program  events 

DESCRIPTION 

This  file  contains  seismicity  parameters  for  one  geographic  area.  It  is  an 
ASCII  file  consisting  of  records  separated  by  new-lmes  The  data  in  the  file  are 
assumed  to  follow  Fortran  77  free-format  rules  (i.e  .  ‘read’  )  Thus,  the  spacing 
between  items  is  not  important,  as  long  as  certain  items,  as  indicated  below, 
begin  a  new  record 

To  assist  the  user  with  documentation,  this  file  may  contain  comment 
fields."  which  are  stripped  out  by  NICE  before  the  Fortran  free-format  reao  is 
called  A  comment  field  begins  with  the  sharp  (*)  character  and  ends  with  end- 
of-record  (new-line  or  end-of-file) 

The  contents  of  the  seismicity  parameter  file  are  as  follows 
Beginnmg-of-file 

*  a-factor  •  ■  b-factor  ■  •  mm-moment  • 

New-line 

-  num-of-depth-pomts  • 

New-lme 


New-line 
<  box -I  at 


box-ion. 


box-iat. 


box-ion. 


<  box-width  > 

<  strike-mean 

<  dip-mean  > 

<  rake-mean 
End-of-file 


strike-sigma  ■ 
dip-sigma  • 

<  rake-sigma  -• 


"New-line"  indicates  where  new-lme  (carriage-return,  characters  are 
required  in  the  input.  Specifically,  the  following  oata  items  must  begm  a  new 


me:  •- num-of-depth-pomts z1  -  and  •  box-iat„  •  For  considerations  cf 


readability,  the  user  is  also  allowed  to  insert  new-nnes  between  other  data  items 
as  well,  as  long  as  none  of  the  three  mentioned  above  is  omitted  As  ment.oned 
above,  commend  fields  are  also  allowed 


The  variables  appearing  above  are  defined  as  follows  Refer 
Appendix  A  for  mathematical  details. 


a-factor 


b-factor mm-moment 


These  are  parameters  of  the  distribution  of  seismicity  with  resoect  to 
seismic  moment  Let  N(M)  denote  the  expected  number  of  eartncuaKes 
occurring  in  the  geographic  area  of  interest,  having  a  seismic  moment  greater 
tnan  M  Then,  this  function  is  given  by 


ogiO  N 


a-factor 


o-factor  •  *  M.  M 
N  -  C  M 


mm-moment 

mm-moment 


V’/. 


•  ■  V  1 
% 


ssa 


■  num-of-deoth-ppints 

This  variable  is  the  number  o'  pomts  cefmmg  tne  ceom  O'Stndut  cn  c* 
seismicity  ■  num-of-depth-pomts  -  must  be  at  least  2  Out  net  create-  mam  25 
[Note  the  maximum  currently  25.  is  given  bv  Fortran  oa-ameter  maxz  ' 
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<z1>  <F1>  <z2>  <F2>  ...  <zn>  <Fn>. 

This  list  is  a  sampling  of  the  distribution  of  seismicity  with  depth.  Tne 
number  of  sample  points,  shown  here  as  'n'  is  identically  <  num-of-deptn- 
points  > .  The  pairs  ( <  z-  > .  <  F(  > ),  i  =  1 .....  n,  are  assumed  to  be  samples 
of  the  function  F(z),  defined  to  be  the  fraction  of  earthquakes  which  occur  below 
depth  z.  By  definition,  the  following  constraints  on  these  values  must  be 
obeyed: 


'  ^  —  u 

<F1  >  =  F(0)  =  1 
<Fn>  =  0 

<  z(  >  >  <  z(  1  >  .  i  =  2,3 . n 

<  Fj  >  <  <  Fj  1  >  ,  i  =  2,  3 . n 

F(z)  is  assumed  to  be  linearly  interpolated  between  sample  points. 

<  box-lat1  >  <  box-ion^,  <box-lat2>.  <  box-lon2  >  ,  <  box-width  >  . 

These  parameters  define  the  spatial  distribution  of  seismicity  for  this 
geographic  area.  Earthquake  epicenters  are  assumed  to  be  uniformly 
distributed  with  a  box-like  region  on  the  earth's  surface. 

EXAMPLE 

Here  is  an  example  of  a  seismicity-parameter  file 

#  Seismicity  parameter  file  for  Hmdu-Kush.  subarea  A 

# 

*  Seismicity  versus  moment  parameters 
» 

1  593  0.565  4,el4  *  a-factor.  b-factor,  min-moment 

tt 

*  Seismicity  versus  depth  parameters 

p 

2  p  number  of  points  in  distribution 

p  z  F 

0  1 
50  0 
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4  Epicenter  box: 

4 

#  latl  ionl  Iat2  lon2  width 

37.4  69.9  39  2  72.0  1  2 

4 

#  Fault  plane  distribution: 

4 


28. 

10. 

4  strike  mean  and  std  dev 

34. 

10. 

4  dip  mean  and  std  dev 

-80. 

10. 

4  rake  mean  and  std  dev 

The  latitude-longitude  pair  <box-lat1>.  <  box-laU  >  are  the  coordinates  ct 
point  shown  in  Figure  A.l.  and  similarly  for  Indices  of  the  points  on  toe 
corners  of  the  box  increase  counter-clockwise  from  P1  The  arc-iengtn  •  bcx- 
width  >  .  the  quantity  2_14  on  Figure  A.l.  is  the  length  of  the  arc  connecting  p. 
and  P4.  The  units  of  these  quantities  is  degrees 

<  strike-mean  >  ,  <  strike-sigma  > 

These  variables  are  the  mean  and  standard  deviation  of  a  uniform 
distribution  for  strike,  in  degrees 

<  dip-mean  >  ,  <  dip-sigma  > 

These  variables  are  the  mean  and  standard  deviation  of  a  uniform 
distribution  for  dip,  in  degrees 


rake-mean  .  <  rake-siama 


FILE  TYPE 


Station  File  .  .  table  of  network  stations  and  their  locations. 


CREATED  BY 
User 


USER  BY 

Program  synwav.  anlyst 


DESCRIPTION 

This  file  contains  the  names  and  geographic  locations  of  the  stations  in  a 
seismic  network  It  is  an  ASCII  file  consisting  of  records  separated  by  new-lines. 
The  file  contains  one  record  for  each  station.  Each  record  contains  four  data 
items  in  Fortran  77  "free"  format  (i.e. .  read  *').  These  items  are  (from  left  to 
right) 

1  Station  identifier  (a  character  string). 

2  Station  code  (an  integer)  (This  code  is  for  future  use.  It  is 
presently  read  and  ignored  It  must  be  included  in  the  file, 
however,  so  as  to  pacify  the  Fortran  77  free  format  utility  ) 

3  Station  latitude  (degrees  north  of  equator). 

4  Station  longitude  (degrees  east  of  Greenwich). 


Additional  items  on  any  station  record  are  overlooked  (and  so  may  be  used  for 
documentation) 
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EXAMPLE 

Here's  a  station  file  for  a  network  of  stations  which  are  either  in  Ethopia 
or  near  large  American  universities: 


AAE  -1  9.0292 

AAM  -1  42.2997 

SCP  -1  40.7950 


38.7656  Addis  Ababa,  Ethopia 
-83.6561  Ann  Arbor.  Michigan 
-77.8650  State  College.  Pennsylvania 
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FILE  TYPE 

Structure  File  ...  a  layered  earth  structure  model  used  for 

surface  wave  calculations 


CREATED  BY 
User. 


USED  BY 

Program  synwav 


DESCRIPTION 

This  file  contains  the  parameters  of  a  plane  layered  earth  model  which  is 
used  by  synwav  for  calculating  surface  wave  responses.  It  is  an  ASCII  file 
consisting  of  records  separated  by  new-lines.  The  data  are  in  Fortran  77  "free" 
format  (i.e. ,  read  *').  Therefore,  arbitrary  white-space  may  be  used  to  separate 
the  items  on  each  record  (including  extra  new-lines) 

The  contents  of  the  file,  record  by  record,  are  as  follows: 

Record  i : 

<  comment  > 

Record  2: 

■  num-of-layers  >  <  umts-flag  > 

Record  3 

<  thickness  1  >  <  a. (  >  <  /?1  >  •-  p 1  >  q/?1  ■  •  dummy  • 

Record  4 

<  thicknesSj  ■  r  a2  '  ‘  @2  ’  '  ^^2  ’  ^ummy  ' 
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Record  2-n: 

<  thickness,,  >  <a  >  <  B  >  <o  >  <  >  <  dummy 

n  n  rn  rr\  ^rn  7 

In  the  above,  subscript  n  stands  for  the  variable  <  num-of-layers  > 
The  meanings  of  the  above  variables  follow. 

<  comment  > 

A  character  string  of  80  characters  or  less 

<  num-of-layers  > . 

Number  of  layers  in  the  earth  structure  model 


<  umts-flaq  > . 

If  <  units-flag  >  is  a  non-zero  integer,  then  the  input  values  of 
<•  thickness!  >  .  <  a  >  ,  <  P{  >  and  <  >  ,  i  =  1 .  .  .  .  n.  are  multiplied  by  1 000  If 

<  units-flag  >  is  zero,  no  action  is  taken.  It  is  an  error  to  have  <  units-flag  >  be  a 
non-integer 

<  thickness(  • 

Thickness  of  the  ith  layer  (meters) 


P 


wave  velocity  of  the  i 


th 


layer  (meters  second) 


S 


wave  velocity  of  the  i 


th 


layer  (meters  second) 


P ,  ' 


Density  of  the  ith  layer  (kg.m3) 
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Shear  Q  of  the  ith  layer 
<  dummy  -> 

A  number  which  is  read  and  ignored  It  must  be  included  so  as  to  pacify 
the  Fortran  free  format  utility 

EXAMPLE 

Here  is  an  example  of  a  structure  file; 

Standard  East  Kazakh  Structure 
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FILE  TYPE 

Table  File  .  tabulation  of  travel-time  and  amplitude  versus 

distance  and  focal  depth,  for  a  particular  body 
wave. 

CREATED  BY 
User. 

USED  BY 

Program  synwav.  anlyst 


DESCRIPTION 

This  is  a  binary  file  containing  variable  length  records.  It  contains  a 
tabulation  of  travel-time  and  amplitude  for  a  given  body  wave  (e  g  .  P)  sampled 
on  a  rectangular  grid  in  focal  depth/epicentral  distance  space 

The  content  of  the  file,  record  by  record,  is  as  follows. 

Record  1 

<  num-of-depth-samples  > 

Record  2 

Record  3 

•- num-of-distance-samples  ■ 

Record  4 

'  *1  *  '  x2  '  xn  ’ 

Record  5 

*  T11  ‘  A1 1  ‘  T21  ’  *  A21  '  Tn1  An1 
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Record  6: 


<T12>  <Ai2>  <T22>  <a22 


Record  4  -  m: 


Tn2  >  <  An2 


<^1m>  <^1m>  <^2m>  <^2m>  '  ^ 

End-of-file 

In  the  above,  subscripts  m  and  n  stand  for  the  variables  <  num-of-depth- 
samples>  and  <  num-of-distance-samples  >.  respectively. 

The  meanings  of  the  above  variables  follow 


<  num-of-depth-samples  > . 

This  is  the  number  of  sample  points  with  respect  to  event  focal  depth. 
(Fortran  type  'integer  .) 


The  value  of  the  r  sampling  depth  (kilometers,  increasing  downward). 
It  is  assumed  that  <  z1  >  is  the  earth  s  surface  and  that  <  z(1  >  <  <  z(  >  . 
(Fortran  type  real  ) 


<  num-of-distance-samples  > . 

This  is  the  number  of  sample  points  with  respect  to  event-station 
epicentra I  distance  (Fortran  type  'integer' ) 


The  value  of  the  r  sampling  distance  (degrees)  It  is  assumed  that 
<x  i  >  <  <x  >  (Fortran  type  real  .) 


I»1 
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<T  > 

_ !L_ 

The  travel  time  (seconds)  at  the  point  <  x|  >  .  <  z^  >  .  (Fortran  type 

real'  ) 

<  A  > 

_ 11 _ 

The  peak-to-peak  amplitude  (microns)  at  the  point  <  > .  <  z  > .  The 

amplitude  table  is  not  required  to  be  normalized  for  an  event  of  a  particular  size 
(see  manual  entry  "Reference  File").  (Fortran  type  real'.) 

Notes 

The  travel  time  and  amplitude  values  at  the  sample  points  are  normally 
non-negative  numbers  A  negative  value  is  interpreted  to  mean  that  the  quantity 
is  undefined. 

The  programs  which  use  this  file  assume  that  the  travel  time  and 
amplitude  at  an  arbitrary  distance  and  depth  are  obtained  through  interpolation 
between  the  given  sample  points.  Values  of  travel  time  and  amplitude  outside 
the  sample  grid  are  not  defined,  i.e..  for  z  <  <  z1  >  ,  z  >  <  z  > .  x  <  <  x1  >  or 
x  >  <  xn  >  Values  are  also  not  defined  in  an  interval  which  is  adjacent  to  an 
"undefined"  (negative)  sample. 

A  discontinuity  in  travel-time  or  amplitude  with  respect  to  distance  or 
depth  may  be  represented  through  the  device  of  zero  width  intervals  (e  g  . 
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APPENDIX  A 

DISTRIBUTION  OF  EARTHQUAKE  HYPOCENTERS. 
RECURRENCE  INTERVALS  AND  MOMENTS 


In  this  appendix,  we  describe  m  detail  the  distributions  usee  to  caicjiate 
source  parameters  for  the  earthquake  catalog  generated  by  toe  EVENTS 
program  These  parameters  are: 

1  Epicenter 

2  Depth 

3  Slip  vector  orientation 

4  Scalar  moment 

5  Origin  time 

A  1  DISTRIBUTION  OF  EARTHQUAKE  EPICENTERS 

Epicenters  are  distributed  over  boxes  on  the  surface  of  the  globe 
Denote  by  =  (0,  0).j  =  i  ,2,3.4,  the  colatitude-longitude  pairs  at  the  corners 
of  the  box,  as  shown  m  Figure  A  1  The  lines  connecting  the  comers  are  great 
circles  and  mtersect  at  right  angles  The  mdeces  increase  counter-clockwise 
from  P.  The  arcs  connecting  P  to  P  have  lenqth  l  and  the  azimuth 

1  3  1  1  3  ij 

(clockwise  from  north)  of  P  from  P  is  a  The  positions  P,  and  P~  and  arc 

I  i  ij  12 

length  ^  2  are  specified  by  the  user 

The  program  initially  tmds  a,2  and  ,rorr  a  stancaro  c  stance- 
azimuth  routine  Then,  if  RV  is  a  random  vanao  e  uniformly  c  stno^.tec  o^e'  jC.il 
the  Quantities 
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are  calculated  A  different  real izat  or  of  the  RV  s  jsed  for  p^  ar'c  c  -  u  T-e 
required  coordinate  pair  PR  =  (0R  0R)  s  given  by 

0R  =  cos'1  [cos  (0r)  cos  (AiR) 

-  Sin  (0^  sin  (A  1  R)  cos  (al2  -  alRl] 

<pR  =  ^  -  sm'  1  [sm  (A  1 R)  sm  !a*2  -  a,R)  sn  (0R)’, 

where 

A.p  =  cos'1  [cos  (pl2)  cos  (p14;J 
fl^P  -  Sin  [ s i n  sin  f  A  ^  R ,'  ] 


A  2  DISTRIBUTION  OF  EARTHQUAKE  DEPTHS 

In  the  program,  the  cumulative  number  of  events  occurring  a:  c'  ce'cw  a 
particular  cepth  is  piecewise  linear  Denote  by  the  va'ues  o4  the  cuTui at  ve 
d ' st r . tout i on  at  depths  Zr  with  Z,  -  0  ispecified  by  the  user)  arc  by  RV  a 
random  var.ab'e  unTormiy  distributed  on  [O.i  j  Toen  the  source  deptn  Zc  ;s 
4Ound  from 


z 


s 


z 


N  1 

I ' 


fZ  -  Z,1 

N  -  N. 


'ere  N 


R 


RV  •  N,  a-'d  ,  sat  s‘:es  N  ■_  NR  ■  N 


A  - 


•VH 

•  v 
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APPENDIX  B 


SYNWAV  FORMULATIONS  AND  NUMERICAL  METHODS 


This  ap perP  x  Tescr  Des  tre  '^'■a:‘'erT‘avca  a*  c^s  a-p 

^et^cas  psep  oy  r^e  SYN/iAV  p'og'ah"  ‘.c  cpr"pJ:e  sy  :*~e?  c  se  s-" 
>s  P  vipea  into  twp  par’s  Pe^ptep  *c  copy  a-'>p  s^ase  Aa<*rs  E  ‘-a;’' 
ea't^cua-e  <s  rr‘CPe  eg  as  a  ppi^r  pp  ,c  e  pc^c  <=■  s:  e  a  r 
h'S!Dry 


B  *  SYNTHETIC  BODv  .VAVE  SFiSMOCRAMS 


T  iwj 


s  the  'ece-ve-  regior'  trans‘er  fjnct'cr  wh.ch  p'ooagates  a  wave 
mQinging  on  the  -eceiver  region  torm  be.ow  to  toe  ea'f~  s  s-dace 
as  eT-ier  vert  ca1  'ad. a-  or  !raosverse  d  sp;acer'e~: 

ft (m,  the  ane  ast>c  attenuation  ope'ator 

!  id]  the  reco'd'hg  instrumentation  'esponse 


For  toe  SYNWAV  module  we  nave  assumed  that  tne  mart  e  ha^s'e' 
,u-  ct  on  obeys  gecmetnc  ray  theory  that  s 


(B"  2; 


//•  ere  R  s  a  gee metr  o  spread  ng  factnr  and  t  is  tne  have  time  The  code 

♦  ^  1 
1  •  os  me  „  i  .,es  c‘  R  an d  t  from  bilinear  interpo  aton  of  tab  es  of  R  ('  z  anp 

G 

1  •  i  ■  s  an  do  ar  distance  arid  z  is  tne  foca1  depth  Tne  .n  te'pciat  on 

./  ■  i-s  ci i s  o  ’  et ur r  t*-e  phase  ve'ccity 
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for  the  P-SV  problem  ana 


S  (w.  6  ^s.  X.  <5) 


for  the  SH  problem,  are  the  displacement  spectra  of  the  up  and  down  going 
waves  at  the  source  The  subscripts  U  and  '  D"  indicate  up  and  dcwngomg 
waves,  respectively  The  orientation  of  the  double-couple  is  specified  by  its 
strike  p  dip  <5  and  rake  X 

We  use  the  formulation  given  by  Aki  and  Richards  ( t 980)  for  the  source 
strength  of  a  point  double-couple  force: 

aP  -  - -  3  [l  •  @  *  ?) 

4^s°s 

aSV  - 1 —  (g  •  1^1  •  2)  (Bi  4) 

4,^s 

a3”  _  — 1— 

4f^s 


wnere  the  vectors  (7  p  f)  form  a  coordinate  system  moving  w  tn  the  array  with 
7  pointing  aiOng  the  ray  and  l  horizontally  and  nc’T.a'  tc  the  d  recton  o‘  ray 
nave  The  tensor  M  >s  the  t  me  rate  of  change  of  the  moment  terser  !•'  a 
coc/d'rate  system  (x  y  7)  whose  un.t  vecto-s  oomt  norm  east  and  ac/m  me 
r riy  vectc'S  a-e 


B-  i 


■V-J-A-V 


.•VV <• 


*  **V>^>*  N  * 


«  . »  » » 


7  =  sin  6  cos  0  x  -  sm  9  sm  0  y  -  cos  6  z 


>< 


p  =  cos  9  cos  0  x  -  cos  6  sin  0  y  -  sin  9  z 
0  =  sin  0  x  -  cos  0  y. 

The  values  of  9  are  computed  from 
6  =  sin'1  (a  /c) 

for  the  downgoing  P  wave  (ag  is  a  at  the  source  and  c  is  the  horizontal  phase 
velocity),  and 

6  =  i-  sm  (a  /c) 

for  the  upgomg  P  wave  Similar  expressions  hold  for  the  S  waves  The 
components  of  the  moment  tensor  in  this  coordinate  are 

M  =  M  (sin  7  cos  X  sin  20  -  sin  26  sin  X  sin  20  ) 

AA  ^  p 

M  =  M  =  M  (sin  6  cos  X  cos  20  -  1/2  sin  26  sin  X  sin  20J 

Ay  yA  u  o  b 

M  =  M  =  -M  (cos  6  cos  X  cos  0  -  cos  26  sm  X  sm  0  ) 

At  U  b  b 

p 

Myy  =  Mo  ^Sin  5  cos  x  s,n  2*s  '  s,n  26  s'n  X  cos  ^s* 

Myz  =  Mzy  =  -  Mq  (cos  6  cos  X  sm  0&  -  cos  26  sin  X  cos  0s; 

Mzz  =  sm  26  sm  X 

where  MQ  -  Mq  m(w).  and  MQ  is  the  scalar  moment  We  have  chosen  a  simq  e 
form  of  the  earthquake  spectrum  m(w)  which  incorporates  the  ejects  on  me  ‘ar 
field  time  history  of  fault  dimension  and  stress  drop  the  fjnct  ora  form  .$ 
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m(u)  = 


(Sip  e- ' 


?) 


IB 


1  ^  ' 


wi  th 


m 


3  »  T 


r.>l 
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The  characteristic  time  r  is  given  by 


r  =  2Ufl 


which  depends  on  the  source  dimension  L  through  the  relation 


L  = 


1/3 


.s. 


5 


where  At  is  the  dynamic  stress  drop  set  to  10  Pascals  MOO  bars)  in  our 
calculations 

SYNWAV  currently  assumes  that  the  source  and  receiver  regon  tr3~s*er 

/  p  \  D 

function  y  (u)  and  T  («)  are  simply  equal  to  one  The  source  anc  rece  ver 
regions  do  not  mnoduce  individual  distortions  to  the  waveform  The  decay  'ate 
(R’  ;  tables  provide  relative  amplitudes  o<  the  phases  P  cP  ana  S 

The  aneiastic  attenuation  operator  is  the  one  per  veo  by  F^tte'^a" 
M962).  given  py 


$ 

$ 


& 


Alw)  e*p  !  ft*  l»  •  2  .n  *  3CC  *: 


where  f  s  freouencv 


:1 
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B  2  SYNTHETIC  SURFACE  WAVE  SEISMOGRAMS 

Because  surface  waves  are  dispersed,  it  is  necessary  to  corrpute  a 
complete  time  senes  at  each  value  of  range  At  a  given  range  source  cepth  and 
orientation,  the  surface  wave  displacement  spectrum  is  computed  from 
(Harkrider.  1970) 


u  ( u  r  ^  ) 


X(“  *  hS]  AR(h,) 


where 


S'*' 


R 


R0  Sin  (r 


--  Sfui) 


exp 


2  •  'k(w)  r  ■  7 r  |  (f  ( w )  I  1  u 


2*f  where  <  s  tfeouency 


the  azimuth  measured  ciockwse  Noon 


tne  source  deptn 


X'w  #  n  a  scwfce  e«.  tat  ct 


tne  s  ).i''  e  so*-  to/’' 


a  t  ..  t »• ;  :  ■  g  .  '  : 

n  s;  i  i  'u  t  i  •  ••a’ 


•  *■'  i 


~  v* 


t'  ;  '  a  e  <» 


'*  ».  ■  r  /  ■  *,i  i .  ►  .  •  : 


"  e  •  ait  a  '  ..i* 


(B2  1 


Ap(w)  =  the  structure  excitation  function. 


I(w)  =  the  instrument  response. 


We  assume  that  for  long-period  teleseismic  surface  waves  the  source  spectrum 
may  be  written  m  the  simple  form 

S(w)  =  -  Mq 

where  M  is  seismic  moment  The  source  excitation  functions  are 
o 

X(ui.*)  -  dQ  -  >[-d1  sm  (ws)  -  d2  cos  (*-*s) 

-  d3  sin  2{f  *s)  -  d4  COS  2(#-*s)  (B2  2) 

where  is  the  strike  (measured  clockwise  from  North)  The  coefficients  d  are 
dependent  on  dip  X  and  rake  <5  and  are  listed  in  Table  B2  1  The  dip  is 
measured  down  from  the  horizontal  and  the  rake  is  measured  counterclockwise 
from  the  horizontal  in  the  plane  of  the  fault  The  eigenfunctions  referenced  m 
Table  B2  i .  are  those  defined  by  Harknder  (1970) 

The  function  <?'(*)  'S  given  by 

f  -  i  for  vertical  component  Rayleigh  wave  (positive  up) 

ie:«)  for  radial  component  Rayleigh  wave.  where  e  ;s  the  eliipticity. 
/  .  for  Love  Wave  (transve'se  component) 
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